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ABSTRACT 

We present preliminary results of an analysis into the effects of primordial voids on the 
cosmic microwave background (CMB). We show that an inflationary bubble model of 
void formation predicts excess power in the CMB angular power spectrum that peaks 
between 2000 < I < 3000. Therefore, voids that exist on or close to the last scattering 
surface at the epoch of decoupling can contribute significantly to the apparent rise in 
power on these scales recently detected by the Cosmic Background Imager (CBI). 
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1 INTRODUCTION 

, One of the primary goals of modern cosmology is to gain 
an understanding of the formation and evolution of struc- 
ture in the universe. Analyses of redshift surveys such as 
the 2-degree field galaxy redshift survey (2dFGRS, Peacock 

■ et al. 2001) suggest that there are large volumes of rela- 
tively empty space, or voids, in the distribution of galaxies. 
It seems that the universe is made up of a network of voids 
with most galaxies tending to be found in two dimensional 
sheets or filaments that surround these under-dense regions. 

In the hierarchical model of structure formation, gravi- 
tational clustering is responsible for emptying voids of mass 
and galaxies (Peebles 1989). Simulations of the standard 
cold dark matter (CDM) model predict significant clumps 
of matter within voids that are capable of developing into 
observable bound objects (Dekel & Silk 1986; Hoffman et al. 
1992). Peebles gives an in-depth discussion of the contradic- 
tions of this prediction with observation (Peebles 2001). He 
argues that the inability of the CDM models to produce the 
observed voids constitutes a true crisis for these models. Ad- 
ditionally, recently announced deep field observations from 
the Cosmic Background Imager (CBI) (Mason et al. 2002) 
show excess power on small angular scales, £ > 2000, in the 
cosmic microwave background (CMB). 

It may be possible to explain the observations by pos- 
tulating the presence of a void network originating from pri- 
mordial bubbles of true vacuum that nucleated during in- 
flation (La 1991; Liddle & Wands 1991; Turner et al. 1992; 
Occhionero & Amendola 1994; Amendola et al. 1996). In 
this scenario, the first bubbles to nucleate are stretched by 
the remaining inflation to cosmological scales. The largest 
voids may have had insufficient time to thermalize before 
decoupling and may persist to the present day. If voids exist 



at recombination they will leave an imprint on the cosmic 
microwave background. On the other hand, if they formed 
much later, their effect on the CMB will be negligible and 
will not be observed with the current generation of experi- 
ments. 

The effects from primordial voids on the CMB have 
been investigated by a number of authors (Thompson & 
Vishniac 1987; Sato 1985; Martinez- Gonzalez et al. 1990; 
Martinez-Gonzalez & Sanz 1990; Liddle & Wands 1992; 
Panek 1992; Arnau et al. 1993; Meszaros 1994; Fullana et al. 
1996; Meszaros & Molnar 1996; Shi et al. 1996; Baccigalupi 
et al. 1997; Amendola et al. 1998; Baccigalupi et al. 1998). 
The most complete investigation was carried out by Sakai 
et al. (1999) who modeled the effect for a distribution of 
equally sized voids. 

In this Paper, we use a simple inflationary bubble model 
to show that if the voids that we see in galaxy surveys today 
existed at the epoch of decoupling, they would contribute 
significant additional power to the CMB angular power spec- 
trum between 2000 < t < 3000. Unlike previous analyses, we 
develop a general method that allows the creation of maps 
and enables us to consider an arbitrary distribution of void 
sizes. We model a power law distribution of void sizes as 
predicted by inflation (La 1991) and also take into account 
the finite thickness of the last scattering surface which sup- 
presses part of the contribution from small voids. 



2 VOID NETWORKS IN THEORY AND 
OBSERVATION 

2.1 Predictions of the inflationary bubble model 

In the extended inflationary model (La & Steinhardt 1989; 
see Kolb 1991 for a review), true vacuum bubbles nucleate 
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during inflation in first order phase-transitions. This model 
predicts a distribution of bubble sizes greater than a given 
radius r of the form, 



Nb(> r) oc r 



(1) 



Typically, extended inflation is implemented within the 
framework of a Jordan-Brans-Dicke theory (Brans & Dicke 
1961). In this case, the exponent a is directly related to 
the gravitational coupling iv of the scalar field that drives 
inflation, 



3 + 



w + 1/2' 



(2) 



Values of Lu > 3500 are required by solar system experiments 
(Will 2001), although models have been proposed that ei- 
ther suppress or hide the present value of u) (La et al. 1989, 
Holman et al. 1990). The main driving force behind these 
models is that a low a can lead to large effects on the CMB 
if arbitrarily large voids are allowed (see Liddle & Wands 
1991 for a review). The normalization of the bubble size dis- 
tribution also depends on ui as well as on the energy scale 
of inflation. 

Once formed, the bubbles will expand and form a shock 
wave on their boundary with the surrounding matter. After 
inflation ends, matter will start to flow relativistically back 
into the freshly created underdensities. However, cold dark 
matter only travels minimally into the void, since it becomes 
non-relativistic early on (Liddle & Wands 1992). Gravita- 
tional collapse of CDM will begin as normal at equality, fur- 
ther emptying any persisting voids. We expect baryonic mat- 
ter to be pushed much further into the void as it is tightly 
coupled to relativistic photons until the epoch of decoupling 
when it will begin to gravitationally collapse back onto the 
CDM. 



2.2 Void detections in redshift surveys 

A number of different void finder algorithms have been de- 
veloped to detect voids in redshift surveys (Kauffman & 
Fairall 1991; Kauffman & Melott 1992; Ryden 1995; Ry- 
den & Melott 1996; El-Ad & Piran 1997; Aikio & Mahonen 
1998). So far, such algorithms have been used to search 
for voids in the first slice of the Center for Astrophysics 
(CfA) redshift survey (Slezak et al. 1993), the Southern Sky 
Redshift Survey (SSRS) (Pellegrini et al. 1989, El-Ad et al. 
1996), the Infra-Red Astronomical Satellite (IRAS) 1.2 Jy 
survey (El- Ad et al. 1997) , the Las Campanas Redshift Sur- 
vey (LCRS) (Miiller et al. 2000), the Updated Zwicky Cat- 
alogue (UZC) (Hoyle & Vogeley 2002) and the Point Source 
Catalogue redshift (PSCz) survey (Plionis & Basilakos 2001, 
Hoyle & Vogeley 2002). These investigations indicate that 
30-50% of the fractional volume of the universe is in the 
form of voids of underdensity 8p/p < —0.9, in line with the 
inflationary model predictions. These voids range in radius 
from r min = 10 h' 1 Mpc to r max = 20-30 h' 1 Mpc. 



3 THE PHENOMENOLOGICAL VOID 
NETWORK MODEL 

We model the voids seen today as spherical underdensities of 
Sp/p = — 1. Each void is bounded by a thin wall containing 



the matter that is swept up during the void expansion. This 
forms a compensated void. We take the background universe 
to be an Einstein-de Sitter (EdS) cosmology, which is a good 
approximation since the majority of the effect on the CMB 
comes from voids on or close to the last scattering surface 
(LSS) and the universe tends towards an EdS cosmology at 
early times. Maeda & Sato (1983) and Bertschinger (1985) 
use conservation of momentum and energy respectively to 
show that these compensated voids will increase in radius 
r v between the onset of the gravitational collapse of matter 
at equality and the present day such that, 

r v (rj) oc rf , (3) 

with j3 ~ 0.39 and where rj is conformal time. 

In this Paper, we consider a phenomenological primor- 
dial void model that is based on the predictions of extended 
inflation with parameters chosen to be in agreement with 
current redshift survey observations; a full analysis of a 
larger family of models will be presented in a subsequent 
paper. Motivated by the inflationary scenario, we assume 
a power-law distribution of void sizes in the universe to- 
day, as given in equation (jl|). We further assume that the 
mechanism creating the voids imposes an upper cut-off on 
the size distribution. A possible mechanism for this cut-off 
could be that the tunneling probability of inflationary bub- 
bles is modulated through the coupling to another field. We 
can therefore go to the limit of large U), leading to a spec- 
trum of void sizes with a = 3. Apart from avoiding problems 
with well-established local measurements of gravity, this as- 
sumption allows us to match the observed upper limit on 
void sizes from the galaxy redshift surveys. 

The minimal present void size is chosen to agree with 
redshift surveys, r m i n = 10 h~ Mpc. For the maximal ra- 
dius, we choose the average value that is found, r max = 
25 ft" 1 Mpc. This scale can be strongly constrained by CMB 
data - if the maximal size were much larger, the voids would 
add too much additional power (given the observed value of 
F v ) at the wrong scales (I < 2000). On the other hand, 
smaller voids would not be able to produce any significant 
contribution to the CMB power spectrum. The exponent of 
the size distribution a is weakly constrained by both the- 
ory and observations. Varying a changes the position of the 
peak and the overall power. This can partially adjust the 
influence of the other parameters (see figure ^). 

We normalize the distribution by choosing the total 
number of voids so as to fill the required fraction of the 
universe today, F v . Redshift surveys point to F v « 0.4, ie. 
40% of the volume of the universe is in underdense regions. 
The positions of the voids are then assigned randomly, mak- 
ing sure that they do not overlap. In order to speed up this 
process, we consider only a 10° cone. This limits our anal- 
ysis to £ > 100, which is satisfactory for our purpose since 
the main contribution from voids is on much smaller scales. 



4 STEPPING THROUGH THE VOID 
NETWORK 

4.1 Voids between us and the LSS 

We ray trace photon paths from us to the LSS for the 10° 
cone in steps of 1'. Each void in the present day distribu- 
tion that is intersected by the photon path is evolved back 
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ip = , for r > r v . 



(9) 



Figure 1. Cross section of a void (Thompson— Vishniac model). 
The trajectory of a photon is depicted from the LSS to an ob- 
server. The subscripts 1 and 2 denote quantities at the time the 
photon enters the void and at the time it leaves, respectively, a 
is defined as the angle formed between the line-of-sight direction 
and the direction of the void's centre. 5a is defined as the scat- 
tering angle of a photon, d is defined as the comoving distance of 
the void's centre, d^gg is defined as the comoving distance of the 
LSS. 

in time according to (^|) to determine whether the photon 
encounters the void. If a photon intersects a void between 
us and the LSS, we compute the Rees-Sciama (1968) effect 
due to the deviation in the redshift of the photon as it passes 
through the expanding void and the lensing effect due to the 
deviation in its path. Thompson & Vishniac (1987) applied 
double local Lorentz transformations at each void bound- 
ary to obtain the redshift deviation (the RS effect) and the 
scattering angle of a photon (the lensing effect), 



ST 

Ars = — 



= (H2R2f COS 9 2 (3/3 - | COS 2 9i 



-0! + e'j + e' 2 -02 = {H2R2) 2 sm{20 2 ) 



(4) 



(•») 



where H is the Hubble parameter, R = ar v is the proper 
length of the void radius, (3 is given by (Eh and the angles 
Ox, 0[, 0' 2 and 02 are defined by reference to figure |lj. 

This treatment agrees with the complementary ap- 
proach of using the potential approximation for the under- 
density (Martinez-Gonzalez et al. 1990). The potential ap- 
proximation is defined in the EdS background making cer- 
tain assumptions for which one of the Einstein equations 
reduces to the Poisson equation, 

-^■V 2 ^ = ^GpS, (6) 
a 1 

where p is the background density and S = Sp/ p is the den- 
sity fluctuation field. For a spherical void with a thin shell, 
p(r], x) is explicitly written as, 



(7) 



p(r], r) = p(v)0(r - r v (rj)) + pi n (ri,r)0(r v (r)) - r) 
+cr(ry)(5Dirac(r - r v {rj)) , 

where is the Heaviside function, foirac is the Dirac delta 
function, pi n is the energy density inside the void, and a is 
the surface energy density of the shell. Using this model of 
a void and assuming pi n to be homogeneous, (|^) is easily 
integrated to give, 



ip — ^-H 2 a 2 (r 2 — r 2 ) 5 , for r < r v , 



The non-linear growth of the void causes the potential inside 
it to grow with respect to the EdS background. This gives 
rise to the RS effect given by, 



Ars = 



(10) 



-2 1 dx ■ Vip . 

'LSS 

The integration of (|l^) for a void between us at the LSS re- 
sults in equation M) , that is the Thompson- Vishniac result 



4.2 Voids on the LSS 

If a photon intersects a void on the LSS, we use the po- 
tential approximation to calculate the Sachs-Wolfe (1967) 
effect due to the photon originating from within the under- 
density. For an empty void (8 = — 1) that satisfies the po- 
tential approximation (8), the SW effect is given by (Sakai 
et al. 1999), 



= (r 2 



= T2#V (r 2 cos 2 6>2 - X 2 ) 



(11) 



(12) 



(8) 



where X is defined as the distance between the centre of the 
void and the LSS. Equation (jll|) takes a maximal value at 
r = corresponding to the case where a photon originates 
at the void centre. 

We take into account the finite thickness of the LSS, 
which suppresses the SW effect for small voids, by averag- 
ing the contribution from a number of photons originating 
from a LSS of mean redshift 1100 and standard deviation 
in redshift 80. We also calculate the partial RS effect (PRS) 
that arises due to the expansion of the void on the LSS as 
the photon leaves it. Using the potential approximation and 
integrating (jw|) we obtain, 

Aprs = - | pRS = 2 ^3 + ^- 

+ 21 g ( -22-) - 2 r2( " 2 -'? i LSs)2co5 M , 

\')LSS/ '?LSS' 7i / 

where tlss and T2 are the size of the void at j^lss and at the 
time the photon leaves the void (772) respectively and the 
angle 02 is defined by reference to figure hi 

Once the photon has reached the last scattering sur- 
face, we know the variation of its temperature as well as its 
position on the LSS and can create a temperature map (fig- 
ure ^) . We then use a flat sky approximation to obtain the 
Ci spectrum of the anisotropies (White et al. 1999; da Silva 
2002) (see figure ^). This figure is the main result of this Pa- 
per and shows that a void model motivated by theory and 
observations can provide substantial power on scales beyond 
t = 2000. 

We point out that primordial void parameters are still 
poorly constrained by both observation and theory. The bot- 
tom panel of figure kI shows a few further example models. 
For a power-law size distribution (as motivated by the infla- 
tionary scenario), large voids become rarer as a is increased. 
Therefore, since void analyses of redshift surveys only sam- 
ple a fraction of the volume of the universe, there may exist 
voids of larger r max than currently observed. Models with 
high r max tend to predict too much power on scales I ~ 1000. 
However, if we take inflationary models with a > 6, as moti- 
vated by eg. Occhionero & Amendola (1994), then the peak 
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Figure 2. The map of the temperature fluctuations on the surface 
of last scattering from the fiducial void model considered. The 
parameters of this model are given by a = 3, r max = 25 h~ 1 Mpc 
and F v = 0.4. 



moves to larger £ and the total power drops. The filling frac- 
tion mainly adjusts the overall power. 



5 CONCLUSIONS 

The cosmic microwave background is an excellent tool for 
probing the distribution of matter from last scattering until 
today. In the case of voids, the strongest signal stems from 
objects at very high redshifts, especially from those already 
present at decoupling. We discuss in this Paper the imprint 
of a power law distribution of primordial, spherical and com- 
pensated voids, which could for example be generated by a 
phase transition during inflation. 

We show that the signature of such a distribution of 
voids, that is compatible with redshift survey observations, 
contributes additional power on small angular scales. At the 
same time, this scenario solves the void-crisis of the CDM 
model. Experiments such as the CBI are able to directly 
probe small angular scales and constrain void parameters. 
We will present a constraints analysis of a wide range of void 
models in a future paper. In this in-depth analysis, we will 
also investigate the non-Gaussian signal of void models that 
are compatible with CMB observations as well as any effect 
of acoustic waves that primordial voids may propagate to 
the sound horizon (Baccigalupi & Perrotta 1998, Corasaniti 
et al. 2001). 

Other sources are also expected to contribute at high 
I. Probably the strongest of these is the thermal Sunyaev- 
Zel'dovich (SZ) effect. Since the thermal SZ effect is strongly 
frequency-dependent, experiments which work at about 30 
GHz (like CBI) will see a stronger signal than those work- 
ing at higher frequencies (Aghanim et al. 2002). Hence a 
mult i- frequency approach should be able to easily disentan- 
gle the contribution of voids from the SZ effect. Unfortu- 
nately, it seems to be difficult at present to predict the 
precise level of the SZ contribution, since different groups 
are reporting different results (see eg. Springel et al. 2000 
for a compilation). Future multi-frequency, high-resolution 
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Figure 3. Top: The CMB anisotropics produced by the fidu- 
cial void model (solid line) compared to the primary CMB 
anisotropics (dashed-triple-dotted). Also plotted are the sum of 
primary and void contributions (dotted) as well as the fluc- 
tuations induced purely by voids on the last scattering sur- 
face (dashed) and by those between last scattering and today 
(dashed-dotted). We show the "standard" cosmological concor- 
dance model: of course a combined analysis of primary and void- 
induced fluctuations would select a different cosmology for the 
primary contribution. 

Bottom: Example models depicting a range of void contributions 
to the CMB fluctuations. The models plotted 
25 ft -1 Mpc and F v = 0.4 (solid line), a = 3, r max = 40/i _1 
Mpc and F v = 0.4 (dotted), a = 3, r max = 40 ft -1 Mpc and 
F v = 0.2 (dashed) and a = 6, r max = 40 h~ x Mpc and F v = 0.4 
(dashed-dotted) . 



and high signal-to-noise maps should be able to significantly 
constrain the contribution of primordial voids to the high £ 
CMB power spectrum. Additionally, deep galaxy redshift 
surveys and measurements of the distribution of matter in 
the Ly-a forest will able to directly explore the presence of 
voids in the baryonic matter distribution at low redshifts. 



ACKNOWLEDGMENTS 

It is a pleasure to thank Andrew Liddle for many interesting 
and crucial discussions. We also thank Antonio da Silva for 
helpful conversations and acknowledge the use of his an- 
gular power spectrum extraction algorithm. LMG is also 
grateful to Charles Lineweaver and the University of New 
South Wales, where some of this work was carried out, for 
their hospitality. LMG acknowledges support from PPARC. 
MK acknowledges support from the Swiss National Science 
Foundation. 



REFERENCES 



Agh 



anjm N. e.t al.. 2002, preprint archived under 



astro-ph/0203112 . 
Aikio J. and Mahonen P., 



1998, Astrophys. J. 497, 534. 



Amendola L. et al, 1996, Phys. Rev. D 54, 7199. 



© 0000 RAS, MNRAS 000, 000-000 



A possible contribution to CMB anisotropies at high £ from primordial voids 5 



1. 

925. 
., 2001, 



Mon. Not. 



Amendola L., Baccigalupi C. and Occhionero F., 1998, Astrophys. 
J. 492, L5. 

Arnau J.V. et al, 1993, Astrophys. J. 402, 359. 

Baccigalupi C, Amendola L. and Occhionero F., 1997, Mon. Not. 

R. Astron. Soc. 288, 387. 
Baccigalupi C, 1998, Astrophys. J. 496, 615. 
Baccigalupi C, Perrotta F., 2000, Mon. Not. R. Astron. Soc, 

314, 1 

Bertschinger E., 1985, Astrophys. J. Supp. 58 
Brans C. and Dicke C.H., 1961, Phys. Rev. 24. 
Corasaniti P.S., Amendola L. and Occhionero F 

R. Astron. Soc, 323, 677 
Dekel A. and Silk J., 1986, Astrophys. J. 303, 39. 
El-Ad H. et al, 1996, Astrophys. J. Lett. 462, L13. 
El- Ad H., Piran T. and Dacosta L.N., 1997, Mon. Not. R. Astron. 

Soc. 287, 790. 
El-Ad H. and Piran T., 1997, Astrophys. J. 491, 421. 
Fullana M.J., Arnau J.V. and Saez D., 1996, Mon. Not. R. Astron. 

Soc. 280, 1181. 

Hoffman Y., Silk J. and Wyse R.F.G., 1992, Astrophys. J. Lett. 
388, L13. 

Holman R., Kolb E.W. and Wang Y., 1990, Phys. Rev. Lett. 65, 
17. 

Hoyle F. and Vogeley M.S., 2002, Astrophys. J. 566, 641. 
Kauffmann G. and Fairall A. P., 1991, Mon. Not. R. Astron. Soc. 
248, 313. 

Kauffmann G. and Melott A.L., 1992, Astrophys. J. 393, 415. 

Kolb E.W., 1991, in The Birth and Early Evolution of the Uni- 
verse, Proceedings of the 1990 Nobel Symposium. 

La D. and Steinhardt P.J., 1989, Phys. Rev. Lett. 62, 376. 

La D., Steinhardt P.J. and Bertschinger E., 1989, Phys. Lett. B 
231 231. 

La D., 1991, Phys. Lett. B 265, 232. 

Liddle A.R. and Wands D., 1991, Mon. Not. R. Astron. Soc. 253, 
637. 

Liddle A.R. and Wands D., 1992, Phys. Lett. B 276, 18. 
Maeda K. and Sato H., 1983, Progr. Theor. Phys. 70, 772. 
Martinez-Gonzalez E., Sanz J.L. and Silk J., 1990, Astrophys. J. 
Lett 355, 5. 

Martinez-Gonzalez E. and Sanz J.L., 1990, Mon. Not. R. Astron. 

Soc. 247, 473. 
Mason B.S. et al., 2002, istro-ph/0205384 
Meszaros A., 1994, Astrophys. J. 423, 19. 
Meszaros A. and Molnar Z., 1996, Astrophys. J. 470, 49. 
Miiller V. et al, 2000, Mon. Not. R. Astron. Soc. 318, 280. 
Occhionero F. and Amendola L., 1994, Phys. Rev. D 50, 4846. 
Panek M., 1992, Astrophys. J. 388, 225. 
Peacock J. A. et al, 2001, Nature 410, 169. 
Peebles P.J.E., 1989, J. R. Astron. Soc. Canada 83, 363. 
Peebles P.J.E., 2001, Astrophys. J. 557, 495. 

Pellegrini P.S., da Costa L.N. and de Carvalho R.R., 1989, As- 
trophys. J. 339, 595. 

Plionis M. and Basilakos S., 2002, Mon. Not. R. Astron. Soc 
399. 

Ryden B.S., 1995, Astrophys. J. 452, 25. 
Ryden B.S. and Melott A.L., 1996, Astrophys. J. 470, 160. 
Rees M.J. and Sciama D., 1968, Nature 217, 511 (1968). 
Sachs R.K. and A.M. Wolfe A.M., 1967, Astrophys. J. 147 
Sakai N., Sugiyama N. and Yokoyama J., 1999, Astrophys. J. 
1. 

Sato H., 1985, Prog. Theor. Phys. 73, 649. 

Shi X., Widrow L.M. and Dursi L.J., 1996, Mon. Not. R. Astron 
Soc. 281, 565. 

da Silva A., 2002, DPhil Thesis, University of Sussex (in prepa- 
ration). 

Slezak E., de Lapparent V. and Bijaoui A., 1993, Astrophys. J 
409. 517. 



Springel V., White M. and Hernquist L., 2001, Astrophys. J. 549, 

681; ibid. Astrophys. J. 562, 1086. 
Thompson K.L. and Vishniac E.T., 1987, Astrophys. J. 313, 517. 
Turner M.S., Weinberg E.J. and Widrow L.M., 1992, Phys. Rev. 

D 46, 2384. 

White M., Carlstrom J.E., Dragovan M. & Holzapfel W.L., 1999, 

Astrophys. J. 514, 12 
Will CM., 2001, Living Rev. Rel. 4, 4. 



330. 



73. 
510, 



© 0000 RAS, MNRAS 000, 000-000 



